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On Lossy Joint Source-Channel Coding In Energy 
Harvesting Communication Systems 

Meysam Shahrbaf Motlagh, Masoud Badiei Khuzani, Patrick Mitran 


Abstract —We study the problem of lossy joint source-channel 
coding in a single-user energy harvesting communication system 
with causal energy arrivals and the energy storage unit may have 
leakage. In particular, we investigate the achievable distortion 
in the transmission of a single source via an energy harvest¬ 
ing transmitter over a point-to-point channel. We consider an 
adaptive joint source-channel coding system, where the length of 
channel codewords varies based on the available battery charge. 
We first establish a lower bound on the achievable distortion. 
Then, as necessary conditions for local optimality, we obtain two 
coupled equations that determine the mismatch ratio between 
channel symbols and source symbols as well as the transmission 
power, both as functions of the battery charge. As examples 
of continuous and discrete sources, we consider Gaussian and 
binary sources respectively. For the Gaussian case, we obtain a 
closed-form expression for the mismatch factor in terms of the 
LambertW function, and show that an increasing transmission 
power policy results in a decreasing mismatch factor policy and 
vice versa. Finally, we numerically compare the performance of 
the adaptive mismatch factor scheme to the case of a constant 
mismatch factor. 

Index Terms —Joint Source-Channel Coding, Energy Harvest¬ 
ing, Distortion, Resource Allocation 


I. Introduction 

W IRELESS sensor networks (WSN) provide a tool to 
gather and disseminate information and thus play an 
important role in unsupervised control systems. In many cases, 
sensor devices are deployed in large numbers to cover a large 
area. As a result, regular maintenance and battery replacement 
for each individual sensor is impractical, if not impossible. 
Thus, to develop truly autonomous sensor networks that do 
not require regular maintenance, it is essential to supply a 
sustainable energy source to each node. 

Energy harvesting (EH), i.e., supplying energy by harnessing 
ambient energy resources such as solar, wind and thermal 
energy, is a promising state-of-the-art solution that can sig¬ 
nificantly improve sensor lifetime. In a typical EH node, 
the energy required for various sensor tasks is incrementally 
harvested from the environment and stored during the course of 
operation. However, due to the stochastic nature of renewable 
energy resources, sensor energy consumption must be managed 
adaptively to achieve good performance. 
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A. Contributions 

In this paper, we focus on the design of data transmission 
policies in EH sensor devices. Specifically, we consider a 
scenario where a single node continuously senses data from 
a source and wishes to transmit this data over a point-to- 
point channel. The harvested energy is stored in a battery 
that may leak energy at a rate which depends on the available 
battery charge. The communication is carried by a joint source- 
channel coding (JSCC) scheme, where in each communication 
block a source sequence of fixed length is mapped into a 
channel codeword whose length depends only on the available 
battery charge. In other words, the mismatch factor between 
the channel symbols and source symbols is adaptive. 

We summarize the major contributions of this paper as 
follows; 

• We formulate and model continuous-time lossy joint 
source-channel coding in an energy harvesting com¬ 
munication system, where the transmission power and 
bandwidth mismatch factor, i.e., the length of channel 
codewords per source symbol, are dynamically adapted 
to the available battery charge. 

• We establish a lower bound on the average distortion 
and show that in the case of infinite battery capacity and 
no leakage it is asymptotically achieved with a constant 
transmission power and a constant mismatch factor. 

• Using a calculus of variations technique, we find achiev¬ 
able locally optimal transmission power and mismatch 
factor policies that minimize the average distortion at the 
receiver. 

• We develop an interesting structural result on the in¬ 
stantaneous distortion. Namely, as long as the battery 
is not depleted, locally optimal transmission power and 
mismatch factor policies will adaptively adjust with the 
battery charge in such a way that the instantaneous 
distortion is maintained to a constant level. 

• Eor a moderate-size battery, we numerically show that our 
proposed scheme with an adaptively varying mismatch 
factor achieves a smaller average distortion compared to 
a scheme with a constant mismatch factor. 

• Eor Gaussian and binary sources, we numerically find lo¬ 
cally optimal power policies and mismatch factor policies, 
both as functions of battery charge. With different leakage 
rates, i.e., zero leakage rate as well as arbitrary non-zero 
leakage rates, we observe that a good transmission power 
policy and mismatch factor policy increases and decreases 
respectively, as the battery charge increases. 
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B. Related Work 

Among prior works that consider lossy communication with 
EH transmitters/receivers are ||2l-|[8l. In ||2l, the mean squared 
distortion of the estimated source symbols at the receiver is 
minimized, where both online and offline EH scenarios are 
considered and the mismatch factor is always one. In 0, 
the problem of energy allocation for data acquisition and 
transmission in WSNs is studied, and the case of a single sensor 
as well as the case of multiple sensors are both considered. 
A similar problem is studied in HI for the case of a single 
source with finite battery and data buffer. Another interesting 
work is 0, where a perturbation-based Lyapunov technique 
is used to obtain an online energy management scheme for 
source-channel coding of correlated sources. In 0, a Gaussian 
source is transmitted over a flat fading channel and the offline 
minimization of the total distortion over a finite-time horizon 
is considered. Therein, the optimal distortion and transmission 
power are found subject to a delay constraint for reconstruction 
of the source symbols at the receiver. In ||71, the problem of 
uncoded transmission over a fading channel is investigated, 
and an optimal energy allocation scheme to minimize the total 
distortion is established. 

Communication systems with EH transmitters/receivers have 
also been studied extensively in the context of lossless trans¬ 
mission. Eor instance, the offline minimization of the trans¬ 
mission completion time for a single source is considered 
in 0, where the battery capacity is infinite. This problem 
has been also extended to the cases of finite battery m, 
multiple access channels im, broadcast channels Ea-m 
and fading channels na. Also in Id, the problem of the 
trade-off between the average queueing delay of packets and 
power consumption of a power grid has been studied using a 
2D Markov chain model. 

In another line of work, throughput/sum-throughput max¬ 
imization is considered for point-to-point channels d, in¬ 
terference channels d and relay channels d, d- A 
more realistic battery model with energy leakage is also 
considered in this context ||2T]| - ||2^ . Specifically in ||23]| . 
offline maximization of throughput in a single-user channel 
as well as a broadcast channel has been studied under battery 
imperfection constraints. In particular, the model assumes that 
only a fraction of harvested energy is stored in the battery 
due to charging/discharging inefficiency, whereas a completely 
different type of battery imperfection, i.e., battery leakage over 
time as in ED, is considered in our paper. 

Another similar work in the framework of lossless com¬ 
munication is ll24ll . where authors have studied online maxi¬ 
mization of long-term average sum-throughput in a multiple 
access channel using the same calculus of variations based 
technique. Despite similarity in the approach taken in both 
works, the nature of the problems are quite different. In the 
current work, we are concerned with lossy source-channel 
transmission in EH communication systems. The lossy nature 
of the problem opens a whole new set of questions that were 
not addressed in ll24l . e.g., how can the trade-off between 
increasing and decreasing transmission power be managed so 
that an optimal average distortion is achieved at the receiver? 
Does an adaptive source-channel coding scheme significantly 


improve the distortion? If so, how exactly should the length of 
channel codes per source symbols be changed in terms of the 
battery state? Our focus here is to minimize long-term average 
distortion in a point-to-point channel that exploits an adaptive 
joint source-channel coding scheme. Different from ll24ll . we 
consider a more general battery model with a non-constant 
leakage rate and have numerically investigated the impact 
of different battery leakages on the compression-transmission 
performance. Moreover, in addition to the transmission power, 
the bandwidth mismatch factor between source symbols and 
channel codewords is also adapted to the battery charge. 

The rest of this paper is organized as follows. In Section HU 
we study the communication model as well as the EH model. In 
SectionUni we present the problem formulation. A lower bound 
on the average distortion as well as the achievable scheme of 
the communication resources are studied in Section |IVl Some 
structural results for the case of Gaussian source are established 
in Section [V] In Section IVIl we provide numerical results and 
simulations. Einally, Section Ivnl concludes this paper. 


H. Preliminaries 


A. Communication Model 


We consider the lossy source-channel transmission of a sta¬ 
tionary memoryless source with general alphabet (continuous 
or discrete) over a point-to-point channel. We first assume 
that the communication is carried over K consecutive blocks 
of joint source-channel coding (JSCC), where the number of 
source symbols in each block is fixed, whereas the length 
of the channel codewords varies from one block to another 
based on the available charge in the battery. Specifically, 
during the z*** block, i = l,...,Ar, a sequence of m inde¬ 
pendent and identically distributed (i.i.d.) realizations of the 
source = (si[l], Sj[2], ■■■, Si[m]), are mapped to a channel 
codeword of length rii, i.e., = {xi[l],Xi\p\, ...,Xi[ni]). 

The input a;"* induces a distribution on the channel output 
vT = 2/iK]), according to the law Py (y”') = 

rifeliIPF| 2 f( 2 /i[fc]|a;i[A:])PA: (a;»[A:]), where '¥Y\x{y\x) is the 
conditional distribution of the stationary and memoryless chan¬ 
nel. At the receiver, after each observation i = 
an estimate s™(y"*) = (si[l], Si[2], Si[m.]) of s™ is made. 

We assume that t is the time duration needed for one symbol 
to be generated/transmitted. Associated with r, Afg = m x r 
is the time duration to generate m source symbols, and /Sia = 
rii X T is the time duration to transmit a channel codeword of 
length rii during the block (see Eig. [D- 


Definition 1. We define the bandwidth mismatch factor be¬ 
tween the z*^ channel codeword block, with duration 
and the source symbol block, with duration Afg, as nita) '■= 
AtcJAts = rii/mfor i = 1, ...,K, where := 
is the channel output time (at the end of the z*^ block) 
that we take as the reference time throughout the paper. 
Likewise, {p{tci)}^i the average power constraints on 
the codewords, i.e.. 


, ni 

— ^\xi[j\\^<p{tci), i = l,...,K. (1) 
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Fig. 1: Consecutive blocks of JSCC, where m and n^, i = 1,..., K, are the number of source symbols and length of channel codewords in 
each block, respectively. 


As demonstrated in Fig. [T] the mismatch factor «:(fci) is 
fixed throughout each block, nevertheless it can change adap¬ 
tively from one block to another. Similarly, the transmission 
power p{tci) can change from block to block. The rate- 
distortion function (i.e., lossy source coding rate) in the 
block, for a source S is given by 

Rs{D,) = mi IiS-,S), ^=l,...,iT, (2) 


where the infimum is taken over all conditional distributions 
Pg|g such that E[d(5, S)] < Di in the block, and I{S\ S) is 
the mutual information between the estimated symbols S and 
the source symbols S. We assume that for a given source S and 
distortion measure d{s,s) the rate-distortion function Rs{D) 
has two thresholds Z?max and Rf^ (where Umax is always finite 
whereas Rf' could be finite or infinite), such that: 

[51] Rs{D) is zero for D > Dmax, 

[52] Rs{D) is strictly decreasing, convex and twice con¬ 
tinuously differentiable for 0 < D < I?max, i e., 
R'siD) < 0 and Rs{D) > 0, and continuous at 
R — Flmax5 

[53] lim^j^o= ~oo , and R's{D) is finite every¬ 
where else, 

[54] lim^ij^o Rs{D) = and Rs{D) is right-continuous 
at £) = 0 if is finite. 

These conditions abstract the form of a rate-distortion function. 
For example, the Gaussian source A/^(0, cr^), with squared-error 
distortion measure d{s,s) = |s — sp has the rate-distortion 
function given by 


Rs{D) = 



0<D <a^ 
D > cr2 


(3) 


and satisfies [S1]-[S4] with Umax = and Rf' = oo. An¬ 
other example is a Bernoulli(p) binary source with Hamming 
distortion d(s, s) = 0 for s = s and d{s, s) = 1 otherwise. The 
rate-distortion function is then given by ||25l, 

H{p) - H{D) 0<D< min{p, 1 - p} 

0 D > min{p, 1 — p}, 

(4) 

where H{p) is the entropy of the source and 


Rs{D) = 


H(D) :=-Dlog^D-{l-D)\og^{l-D). 


(5) 


For the binary source, I?niax = min{p, 1 —p} and Rf^ = H{p). 

In the block, the transmit power is limited to pita), and 
thus the channel coding rate i?c(p(fci)) given by 

Rcipitci)) = sup I{X;Y), i = (6) 

where the supremum is taken over all channel input distribu¬ 
tions Fx that satisfy the power constraint p(tci) in ([I]l for the 
block. This model allows for fast fading, where the block 
length is sufficiently large such that the effect of fading can be 
averaged. We assume that the channel coding rate Rc{p) has 
the following properties: 

[Cl] Rc{p) is strictly positive for p > 0, zero and right- 
continuous at p = 0, 

[C2] Rc{p) is strictly increasing, concave and twice con¬ 
tinuously differentiable for p > 0, i.e., R'^{p) > 0 
and i?"(p) < 0. 

Similarly, these conditions abstract the form of a channel 
coding rate function, and an example is the Shannon rate 

function Rc{p) = -^^og 2 {l +p/N), which will be assumed 
later in the paper. 

We assume that the block lengths m and rii, i = 1,..., AT, 
are sufficiently large that the rate-distortion function and the 
channel coding rate have operational significance. Since the 
mismatch factor is fixed throughout each JSCC block, source- 
channel separation holds in each block. Therefore, based on 
a separate source and channel coding scheme, the relation 
between source coding rate Rg{Di) and the channel coding 
rate Rc{p{tci)) is given by 

Rs{Di) = K{tc,)Rc{p{tci)), i = (7) 

One should note that by the condition [S2] and continuity of 
Rs{Di) at Di = Umax, the inverse function Di{Rs) of the 
rate-distortion always exists for 0 < K{t(,i)Rc{p{tci)) < R^^ 
(though in most cases there is no closed-form), and we can 
compute the distortion Di{Rs) in the block in terms of 
the transmission power pita) and the mismatch factor K^ta) 
using (|7]i, i.e., Di{Rs) := U(p(fci), Moreover, if Rf^ 

is finite then the distortion is D{p{tci), = 0 for pairs 

p(tci),K{tci) such that Hi{tci)Rc{p{tci)) > A**'. 
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From the definition of the mismatch factor and Fig. [T] we 
directly obtain 

^ 1 

Similarly, the total average distortion per source symbol over 
K blocks can be written as 


K 


KAt. 


■ ^ D{p{tci), K{tc,))Ats 


(9) 


2=1 

K 


2=1 ^ 


B. Continuous-Time Model 

As a practical assumption, we suppose that both the source 
and channel block times, Atg and Ata respectively, are 
small compared to the battery dynamics. Hence, we develop a 
continuous-time model in the asymptotic regime based on ([Sll 
and (fTOl) . and let K grow so that 


Ts := KAtg 
K 

Tc := ^ Atci , 
i=l 


where the terms Tg and Tc have fixed values as Atg 0, 
Atci —?> 0 and K oo. Eq. ® has the form of a Riemann 
sum in the limit of Atg 0, At^ 0 and K ^ oo, and 
thus the continuous-time limit takes the following form 


We define 



p{Tc) ■■= 


h 

Tc 




( 11 ) 


( 12 ) 


In this paper, we are interested in infinite-time horizon com¬ 
munication and long-term average distortion. We thus let the 
transmission time become asymptotically large, i.e., Tc oo. 
If limr^-i-oo p{Tc) > 1 or limT„^.oo p{Tc) < 1, the backlog 
in either the source symbols queue or the transmission of 
channel codewords tends to infinity. Thus, for a stable source- 
channel communication system we assume limr^-s-oo p{Tc) = 
limT^-yooTg/Tc = 1. Therefore, from (fT^ for an asymptot¬ 
ically stable joint source-channel communication system we 
require that. 


lim — 

Tc^oo T 


c Jo 


K(tc) 


dtr = 1 . 


Hence, we rewrite the continuous-time limit of (fTOl) 


(13) 


-^avg= lim — / D{p{tc),K{tc)) 


Tc-J-oo T. 


s Jo 


K{tc) 


dtc 


= lim — 
T„->.oo T, 


c Jo 


D{p{tc),K{tc))—^dtc, (14) 

K{tc) 


where (fl4li follows from the fact that limT,,->oo Tg/Tc = 1. 


C. Storage Model 

We assume that the energy arrival times and their amounts 
are not known at the transmitter a priori. Therefore, the 
instantaneous energy of the battery is a stochastic process that 
can be characterized in terms of the energy arrival process 
and battery depletion rate. Let {L'ilti’ denote the size of the 
energy packets arriving to the battery at time instants 
where > 0 and Ti < < ■■■ ■ We assume that energy 

packets are i.i.d. with the tail distribution function denoted by 
B{z) = P[i? > z], and the corresponding arrival times are 
a homogeneous Poisson Point Process (PPP). Specifically, the 
inter-arrival times are i.i.d. and exponentially distributed with 
parameter S, i.e., Ar„ := r„+i — T„ ~ Exp((5). One should 
note that the PPP assumption can subsume cases of EH systems 
that are modeled as bursty packet arrivals at discrete times such 
as regenerative shock absorbers that use the piezoelectric effect 
to transform random mechanical shocks into electrical energy. 
Due to the unpredictability of their occurrence, a PPP is a good 
model to describe the behavior of such EH systems. 

The total harvested energy at the transmitter up to the time 
t, {A(f) : f > 0} is thus a compound Poisson process given 
by 

A(f) := (15) 

nSN 

We assume that the capacity of the battery L is finite, that is, 
L < oo. Eurthermore, we assume that the battery is imperfect 
in the sense that when it is charged, it leaks energy over time at 
a rate which depends only on the current battery charge Z{t), 
and is denoted by i{t) at time t. Also, it is clear that there is 
no leakage when the battery is depleted. 

The instantaneous battery charge at time t, {Z{t) : f > 0} 
is therefore a stochastic process described by 

Z{t) = zo + A(t) - J (p{s)i{s)^ds - R{t), (16) 

where Zo = Z{t)\t=o is the initial battery charge and R{t) is 
a non-decreasing, non-negative and continuous-time reflection 
process with R{t)\t=o = 0, that only increases over the set 
{t : Z{t) = L} ll26ll . The reflection process accounts for the 
excess energy arrivals that overflow the battery capacity and 
ensures that even for large energy packet arrivals the storage 
process is bounded by its capacity limit at all times, i.e., Z{t) € 
[0, L]. Eurthermore, p{t) + £{t) is the instantaneous battery 
depletion rate. Since the transmission power p{t) is adapted 
to the battery charge Z{t), this depends on time only through 
Z{t), i.e., p{t) i{t) = p{Z{t)) -f £(Z(t)). Likewise, the 
instantaneous mismatch factor K.(t) is modulated by the battery 
charge Z{t), where K(t) = K(Z{t)). More specifically, since 
p{t), i{t) and K{t) all depend on t only through the battery 
charge Z{t), with slight abuse of notation we denote by p(z), 
£{z) and k{z) the explicit dependence of these on z. In the rest 
of this paper, we refer to p(z) and k{z) as the power policy and 
mismatch factor policy. The storage process Z{t) can then be 
viewed as a continuous-time Markov process, where the state 
space of the process is the finite interval [0, L]. We impose the 
following conditions on the feasible set of power policies and 
leakage rates 

. Vz e (0, L], p{z) >0, and p{z) \ = 0, 
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• SUPo<;,<i,p(2:) < oo, 

• supo<^<i4z) < oo, and = 0. 

The first condition is to avoid a reserve of energy in the 
battery that can never be consumed by transmission and thus 
effectively reduces the usable energy stored in the battery. The 
second condition reflects the fact that instant depletion of an 
amount of energy AE > 0 is not allowed. Thus, an optimal 
power policy p{z) must satisfy both of these constraints. With 
these conditions on p{z) and £{z), Z{t) becomes irreducible 
in the sense that there is only one single communicating class 
in the state space. 

We define ^tAz) as the empirical distribution function of 
the storage process with respect to the reference time, i.e., 

1 

^{Z{ta)<z}dtc. (17) 

Jc Jo 

By the strong law of large numbers, as Tc —?► oo, 
converges to the stationary probability measure of the storage 
process denoted by Tr{z}, almost surely for every value of z. 

The following theorem identifies stationary distribution of 
the storage process Z{t) lIZTll . where we recall that the inter¬ 
arrival times are i.i.d. and exponentially distributed with pa¬ 
rameter 6, i.e., Tn+i — Tn ^ Exp(i5). 

Theorem 1. For L < oo, the storage process Z{t) is positive 
recurrent and there exists a unique stationary probability 
measure tt{z) = P[Z(f) < z], which may have an atom 
ttq := 7r(z)|^=o > 0> is absolutely continuous on (0, L] 

such that ^ 

tt{z) =tto+ [ f{u)du, (18) 

Jo+ 

where f{z) is the absolutely continuous part of the probability 
measure. Moreover, 

f{z){p{z)+Az)) = 5TToB{z)+5 [ B{z-u)f{u)du. (19) 

Jo+ 

Remark 1. From (117b . ttq is the fraction of time that the battery 
remains discharged. 

Remark 2. Equation (119b is the equilibrium condition between 
the rate of down-crossing f{z){p{z) -|- £(z)) and the rate of 
up-crossing SttoB(z} -h S B{z — u)f{u)du at level z. 


7^0 


f{z) 


dz (23) 


k (0) Jo+ k{z) 

= 1, (24) 


where k(0) := k(z)|^=o is the number of channel uses per 
source symbol when the battery is exhausted. Similarly, we 
apply the ergodicity argument to the distortion function in (fT4b . 
We first rewrite 79avg from (fl4t as 


T^avg — 


1 

lim — / D(p{Z{tA),K{Z{tA)) 

Tc^OoIcJq \ / 


1 


k{Z (fc)) 


dtc- 

(25) 


Due to the ergodicity of the storage process Z{tA, we obtain 


79avg — 


[ D {p{z),k{z)) ^-n{dz) a.s. (26) 

Jo Am 

= 7roi9l(p(0),K(0))-b [ dAp{z),kA)) fiAdz 


= 


Dl (piZ), 


AZ) 


(27) 

(28) 


where 


(p,-\.= D{p,k)-. (29) 

\ K J K 

We now wish to find an optimal transmission power policy p{z) 
and mismatch factor policy AA such that the distortion at the 
receiver is minimized. More specifically, we want to minimize 
(l27l) subject to the constraints in (l20b . (l2Tb and (l24b . i.e.. 


,, , inf TroDAp{0),AA)+ [ DApA), AA) AAdz 

f{z),77o,K{z},K{0} Jq+ 


(30) 


s.t. f(z) {p{z) + £(z)) = (tto + e^^f{u)du) (31) 

770 + / f{z)dz = 1 
Jo+ 


ttq 

k(0) 


/o+ 


fiA 

AA 


dz = 1 


f{z) >0, ttq > 0, AA > AA > 0) 


(32) 


(33) 

(34) 


We now assume that the packets of energy are exponentially 
distributed with parameter A, i.e., B{z) = exp(—Az). There¬ 
fore, ( fT^ reduces to 


fiA ipiA + AA) = de ttq -b 


/o+ 


e^“/(u)dt 


( 20 ) 


III. Problem Formulation 
Since, 7r(z) is a probability measure we have 

pL pL 

/ TT{dz) = TTo -b / fiz)dz = 1. (21) 

Jo Jo+ 

By combining (O and (HI, we also obtain the following 
constraint on the mismatch factor 

rTc 


1 


lim — , , „, ,, 

Tc,—^oo Tc Jq AZAc)) 


dtc = Eti 


1 


.AA. 


a.s. 


( 22 ) 


where the strict inequalities AA > 0 und k( 0) > 0 are implied 
by (l33b . We note that for a given /(z), p(z) can be calculated 
directly using (OTb . 

We note that in the formulation of the optimization problem 
(l30b . the mismatch factor and the transmission power can 
be viewed as two limited resources that must be allocated 
efficiently. Specifically, due to (l33l l. the inverse mismatch factor 
has an average of unity. Consequently, using a large mismatch 
factor, which is desirable as this reduces the instantaneous 
distortion, must be balanced with using a small mismatch factor 
to maintain the average constraint. Similarly, the transmitter 
cannot consume more energy than what is stored in the battery, 
which itself is replenished at a maximum rate of b/A (ignoring 
the energy that is lost due to overflow). Thus, periods with a 
large transmission power must be accompanied by periods with 
a small transmission power. 
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IV. Main Results 
A. Distortion Lower Bound 

In this part, we derive a distortion lower bound using 
Jensen’s inequality and convexity of the distortion function. 

1) Finite Battery Capacity: We first note that for any power 
policy p{z), when the battery capacity is finite, the following 
upper bound has been computed in ll24l : 

E. [p{Z)] < ^ (1 - e-^^) . (35) 

One should note that even with a non-zero leakage rate, the 
upper bound in (l35l l still holds, however it is a potentially 
looser bound. 

We now find a lower bound on the objective function 
E^[U^(p(Z), 1 /k(Z))], based on the upper bound in dlSl l and 
the following convexity lemma. 

Lemma 1. The function D'^{p,q) defined in (129b is jointly 
convex over the pair p and q, where q := 1/k. 


Proof: See Appendix A. ■ 

Based on Lemma [T] as well as Jensen’s inequality, we 
establish a distortion lower bound as below 




p(Z),. 


<Z))\ 


> uM E^ [p{Z)] ,E^ 


<Z)\ 


(36) 



:= -Dlb, 


where (EZli follows from (l35]) and the fact that the function 
dKpM k) is a non-increasing function of p, and further 
recalling that E.n.[l/K(Z)] = 1 from (l24b . The lower bound 
in dJTb holds for any transmission power p{z) and mismatch 
factor k(z) that satisfies dH and dSli. 

2) Infinite Battery Capacity: When the capacity of the 
battery is infinite, i.e., L —?> oo, and therefore the chance of 
battery overflow is zero, for any ergodic power policy p{z) the 
lower bound in (IJTI i simplifies to 


E 


TT 


(^p{Z), 




(38) 


One should note that if the leakage rate is zero, every ergodic 
power policy p[z) is such that E^r [p{Z)] = S/X. This suggests 
that in the infinite battery capacity case with no leakage, 
separate source and channel coding with a constant mismatch 
factor k{z) = 1 is optimal. More precisely, consider the 
following choice of power policy 

p*(z) = -+e, z>0, (39) 


with p*(0) = 0, and e > 0 is a small positive number that 
ensures the storage process is positive recurrent, and choose 
K*{z) = 1 for z > 0. Then, 

+ e)(l-7ro), (40) 


E4p*(^)]= 


and since we must have E^[p*(z)] = S/X, then 


S/X 

S/X + e' 


In particular, we then compute the total average distortion of 
this scheme as 


E„ 


U'l' (^p*{Z), -^7^^ =(l-7ro)i9''' -f e, l^+TToUmax- 

(42) 


As e -A 0 the atom ttq tends to zero from (1411 1 and the lower 
bound in (|3^ is thus asymptotically achieved. 


B. Achievable Resource Allocation Scheme 

The two functions p{z) and k{z) are the policies that we 
wish to design in this subsection to manage the limited re¬ 
sources of the system in such a way that the average distortion 
in is minimized. To this end, we use a calculus of 
variations technique which provides necessary conditions for a 
local and therefore global optimal solution to our optimization 
problem. 

1) Calculus of Variations: We define f^{z) and 1/k'^{z) as 
a perturbed density function and perturbation of the inverse 
mismatch factor, respectively, i.e.. 


:= f{z) + eh{z), 


(43) 

(44) 


k‘^{z) k{z) 

where h{z) and g{z) are continuous and bounded perturbation 
functions on (0,L], with h(0+) = h{L) = 0 and p(0“'') = 
g{L) = 0. For sufficiently small e > 0, the perturbed density 
function f’^{z) satisfies (l32b only if 

p L pL pL 

TTo+ f''{z)dz = ttq + / f{z)dz + e / h{z)dz (45) 

c/o+ j0+ J 0+ 


= 1 , 

which due to (ISTb is true for all e > 0 iff 


f h{z)dz = 0. 

J 0+ 


(46) 


(47) 


In addition, from (l33T l we derive the two following conditions 

Tip 

k(0) 

TTp 


k(0) 


/p+ K,^[Z) 

(48) 

'p+ Hz) 

(49) 


We simplify (l48b using (l44b as follows 


TTp 


-f 


(0) 7o+ \K,{z) 


+ ^9{z) f{z)dz 


TTp 


f{z) 


k{0) Jo+ k{z) 

= 1 . 


dz ■ 


/o+ 


g{z)f{z)dz (50) 
(51) 


Similarly, we simplify ( l49b using ( l43b as 

nz) 


Tip 

k(0) ' Jo+ k[z) 

TTp , 


■dz 


fiz] 


k(0) 

= 1 . 


/p+ 


k{z) 


dz 


/p+ 


Hz) 

Hz) 


dz 


(52) 

(53) 


1 - TTo 


(41) 
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E.. 

= 7roi:)'l'(p(0),K(0)) 


D{p{z),k{z)) 


dD{p{z),K{z)) dp^[z) 


= Ett 


U(Z),. 




(^D{p{z),n{z)) 


dp{z) de 

Hz) dD{p{z),Hz)) dp\z) 

HH dp{z) de 


X1 & 

e=0y V Hz) 


/(^ 


; = 0 Hz) 


dz. 


(63) 

(64) 


Thus, analogous to the constraint in (l47l i. both (fSOl l and 
are true for sufficiently small e iff the perturbation functions 
also satisfy the following constraints 

/ 9 {z)f{z)dz = 0 (54) 


/o+ 


Hz) 


dz = 0. 


(55) 


Jo+ Hz) 

For sufficiently small e, two necessary conditions for a local 
and therefore a global optimal solution to the optimization 
problem in (I30l) -(l34li are 


E, 


i^p^iZ), 

(p{Z),- 


1 


HZ) 

^) 


> Ett 

> E, 


(piZ), 

(p{Z), 


1 


HZ)) _ 

1 


K. 


{z))\ 


, (56) 
, (57) 


where Ett® is the expectation with respect to the probability 
measure with the perturbed density function f^{z), (see (l43l l'). 
We now expand the l.h.s of (l56T l as below 


E^ 


(p^iZ), 


1 


Hz)J\ 

= Tro-Dt (p(0),k(0)) + / D\p^{z),Hz))t{z)dz. (58) 
Jo+ 

We then use OTl i to compute p'^{z) as follows 
P"{z) 

, (tth + 

(59) 


+ t(,) 


fiz) 

^ (^0 + fo+ H^f(u)du + ef^l e^^HHdu) _ 


f(z) + eHz) 


(60) 


We also compute the derivative of p'^(z) with respect to e as 
follows 


dpHz) 


de 


e=0 


f(z) 




( 61 ) 


Based on dMl l. we expand D{p^{z), k{z)) to first order in e, 
i.e., 


D{p%z),Hz)) 

= D{p{z),Hz)) 


dD{p{z), Hz)) dp^{z) 


dp{z) 


de 


e=0 


oiH)- 

(62) 


Substituting D{p'^(z), k(z)) and f^{z) into (l58l l results in (l63t - 
(HHl, where we have neglected the higher order terms of e 
(i.e., 0{e^)). By substituting (l64l i into (l5^ . we establish the 


following necessary condition for a local and therefore a global 
optimal power policy p{z) 


D{p{z),Hz))^^^dz - 


dD(j){z),Hz)) Hz) 


Jo+ 

+ <5 

= 0 . 


Hz) 

dD{p{z),Hz)) 1 

0+ 


0+ 


p{z)dz 


dp{z) Hz) \Jo+ 


dp{z) Hz) 
^-H-^)Hu)du] dz 


(65) 


By changing the order of the double integral in the third term 
of (l65T l we obtain 



i{u) 


dp{u) 

dD{p{z),Hz)) p{z) 
dp{z) Hz) 


-du 


D(j>{z),Hz)) 


i(z) 


Idz = 0. 


( 66 ) 


Equation ( |66] | holds for all perturbation functions h(z), that 
satisfy (EZll and (l55l l. We can thus rewrite (l66l l as follows 




dD{p{u),HH)^ 

dp{u) Hu) 

dD{p{z),Hz)) p{z) , Cl 


-du ' 


dp{z) Hz) Hz) 


+ Cn 


D{p{z),Hz)) 

Hz) 

dz = 0, (67) 


where Ci and C 2 are two free parameters. Therefore, based on 
the fundamental lemma of the calculus of variations we derive 
the following necessary condition 


Se^^Hz) 


dD{p{u), Hu)) e 


dp{u) Hu) 


du + D{p{z),k{z)) 


dD{p{z),Hz)) 

dp(z) 


p{z) + Cl + C 2 k{z) = 0, 


( 68 ) 


where z > 0. Eq. ( l68T l is an integro-differential equation 
involving p{z) and k{z). However, by multiplying both sides 
of ( |68] | by and differentiating both sides with respect to 
z and further simplifications, we derive a first order non-linear 
autonomous ordinary differential equation (ODE) equivalent to 

(|6i. 

We now consider the inequality in (l57l i. where the l.h.s can 
be written as 


E, 


DUp{Z),- 


n^iZ)J\ 

= ttoDHp{0),kH)) + f D'<{p{z),K^{z))f{z)dz. (69) 
J 0+ 
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E, 

= E. 


(p{Z), 


1 


dUp{z), 


K-[Z)J\ 
1 


n{Z)J\ 


D{p{z), K{z))g{z)f{z) - j ^ 2 :. (73) 


/o+ 


dK{z) 


We compute k'^{z) using (l44l) as follows 

K^iz) = -, (70) 


where its derivative with respect to e at e = 0 can easily be 
obtained as below 


—^ = -g{z)n^iz). (71) 

ae e=o 

Similar to (IbZt . we write the Taylor series of D{p{z), k^{z)) 
up to first order term in e as below 


D{p{z),K%z)) 

= D{p{z),k{z)) + 


dD{p{z), k{z)) dK'^{z) 


dK{z) 


de 


e=0 


0(e2). 

(72) 


Therefore, by substituting (l44li . (ItTI i and (f72l i into dfi^ . this is 
then reduced to (|7^ . 

By substituting (l73l l into (l57l i and neglecting the higher order 
terms of e, we establish the following necessary condition for 
a local and therefore a global optimal mismatch factor k,{z). 


/o+ 


9{z)f{z) 


D{p{z),K.{z)) - 


dD(p{z),K.{z)) 
dK{z) ‘ 


■,{z) 


dz = 0. 
(74) 


Equation (|7^ holds for all perturbation functions g{z) that 
satisfy (l54l i. Therefore, we rewrite (f74l i as 


9iz)fiz) 


10 + 


D{p{z),k{z)) - 


dD{p{z),K{z)) 

dK{z) 


iz) + /3 


dz= 0, 
(75) 


where /3 is a constant. Based on the fundamental lemma of the 
calculus of variations, we thus derive the following necessary 
condition 


D{p{z),k{z)) 


dD{p{z),K{z)) 


k{z) + (3 = 0, 


z > 0. 
(76) 


Solutions to the equations in (1681) and (l76l) determine a lo¬ 
cally optimal power policy p{z) as well as a locally optimal 
mismatch factor policy k{z). 

2) Instantaneous distortion: In this part, we discuss an 
interesting consequence of (|7^ on the instantaneous distortion 
D{p{z), k{z))^ z > 0. Specifically, as formally presented in 
Lemma ^ a locally optimal solution has the property that 
the power policy and the mismatch factor policy are adjusted 
in such a way that D{p{z), k{z)) is constant for z > 0. 
In other words, if the transmission power is decreased (or 
increased) due to a change in the battery charge, a locally 
optimal mismatch factor will always dynamically adjust so 
that the instantaneous distortion is maintained to a constant 


level. As elaborated in section mil one can think of the 
transmission power p{z) and the mismatch factor k(z) as 
limited communication resources that must maintain long¬ 
term averages, and thus the transmitter is trading off one 
for the other. More precisely, since the mismatch factor must 
satisfy d24li (i.e., the inverse mismatch factor averages to one), 
as the availability of one communication resource (say the 
transmission power) increases due to a large battery charge, 
the transmitter employs a large transmission power and saves 
on the other communication resource (i.e., the mismatch factor) 
by then using fewer channel uses per source symbol. Likewise, 
when the battery charge is low, the transmitter reduces its 
transmission power, but employs a large mismatch factor. 


Lormally, we have the following Lemma. 


Lemma 2. 

/5 /?max> 

The constant fd in (176b must be in the range /3min < 
where 


Pmax ■— lini T-it f -r\\ ’ 

Dio R'^{D) 

(77) 

and 

Mmin ■ — Imi , , -^max- 

i’ti’max R'siD) 

(78) 


Furthermore, for every such a choice of jd, dZU) results in a 
unique constant solution for D{p{z), k(z)), z > 0. 


Proof: See Appendix IbI ■ 

It is easy to verify that for both Gaussian and binary 

r n ur 

sources, lim£, 4 .o = 0 and 


R'AD) 


R'siD) 


= 0 . 


Therefore, for the Gaussian source we have the bound —tr < 
/? < 0. Likewise, for the binary source we have the bound 
— min{p, 1 — p} < /3 < 0. 


Remark 3. Note that the derivation starting with (I691 l and 
leading to the differential equation for the mismatch factor 
in dzSli has not used the assumption that B{z) = exp(—Az). 
Thus, the structural result on the instantaneous distortion being 
constant for z > 0, Lemma 2 as well as the mismatch factor 
and instantaneous distortion in the Gaussian case to be shown 
in dSB and ( 195b are valid for all energy distributions B(z). 


C. A Constant Bandwidth Mismatch Factor 

So far, we have studied a general JSCC scheme where the 
mismatch factor is adaptively adjusted according to the avail¬ 
able battery charge. However it is also interesting to compare 
the results with the simpler scheme where the mismatch factor 
is fixed. Thus, we now consider the case where the mismatch 
factor is constant and does not adapt to the battery charge. Lor 
a fair comparison with the general case of dynamic mismatch 
factor, we retain the constraint (l3^ which results in a constant 
bandwidth mismatch factor of unity, i.e., k{z) = 1, Vz > 0. 
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Therefore, there is only one design parameter, the transmission 
power p{z), to minimize the total average distortion at the 
receiver. We thus have 


Rs{D{z)) = Rc{p{z)), (79) 


and thereby the distortion-rate function is computed only 
in terms of the transmission power p{z). The optimization 
problem is therefore described as follows 


inf 7roi9(p(0)) + / D {p{z)) f{z)dz (80) 

J0+ 

■t- f{z) {p{z) + ^{z)) = (no + e^'^f{u)du^ (81) 

7^0 + / f{z)dz = 1 (82) 

Jo+ 

f{z) >0, ^0 > 0, (83) 


where D{p{z)) := D{p{z), 1) is the average distortion with 
bandwidth mismatch factor of unity. Moreover, the distortion 
lower bound in this case is the same as Dlb in (l35T l. As a 
necessary condition for a local and thus a global optimal power 
policy p{z), we can directly obtain the following equation by 
replacing k{z) = 1 into (l68l l with the substitution A(C'i + 
C 2 ) = C, 


S f D'{p{u))e du+D{p{z))—D'{p{z))p{z)-\-C = 0, 

J Z 

(84) 


where D'{-) denotes the derivative of D{-), and C is a free 
parameter. As explained below dbSI) . we can derive a first order 
non-linear autonomous ODE equivalent to (l84l i which is found 
to be 


Moreover, the lower bound on the average distortion in this 
case is 


-1 


( 88 ) 

To determine the structure of a locally optimal scheme, we first 
replace D{p{z), k{z)) in (l76l l by its closed form expression in 
dSbl l to obtain 


IE„ 


p(Z), 


k(Z) 


> cr^ 


A ]V 


which can be rewritten as 



+/3 — 0, z > 0, 

(89) 


— 1 — k(z) In 


p(z) 


exp (—1 — k{z) In 




(90) 


for z > 0. Therefore, we have 


Wr,. 


k(^;/3) = -- 


/3 


In 1 + 


Pjz) 

N 


z > 0, (91) 


where Wn{-} denotes the LambertW function 1281 that takes 
either real or complex values and has an infinite number of 
branches, each denoted by an integer n. The notation k (z; /3) 
emphasizes the dependence of k(z) in IMl on the choice of 
/3. Moreover, dMl shows that k(z; (3) is a decreasing function 
of z, whenever p{z) is an increasing function of z and vice 
versa. 


XD(j){z))+{5 - Xp{z)) D'{p{z))+p{z)p’{z)D’'{p{z))+C = 0, 

( 85 ) 

where £>"{■) denotes the second derivative of D{-). 

Remark 4. This ODE is identical inform to (102) in 
Any solution p{z) of the ODE in (l85l for C < —XD{j) is 
non-decreasing in z, for p{z) > 0. The proof of this is similar 
to that of 4241 Lemma 1] except for the change in the direction 
of the inequality which is due to the fact that the rate function 
r(jp) in (102) of 1 1241? . which is concave, is replaced with the 
distortion function D{p) in ( 185b , which is convex. 

V. Gaussian Source and Channel 

In this section, we specialize our results for a Gaussian 
source using the rate-distortion function Rs{D) given in (O. 

We assume that the Shannon rate function Rc{p) = - log 2 (l + 
p/N) is considered for the channel coding rate. Consequently, 
the distortion function D{p, k) can be computed as 

D{p{z),k{z)) = a'^ , ( 86 ) 

or equivalently 

iAt(p(z),lMz)) = ^(^l + ^) . (87) 


Remark 5. Based on the properties of the LambertW func¬ 
tion, n = 0 and n = —1 are the only branches that yield a 
real value for Wn(x), where Wq{x) > —1 for x > —1/e and 
W-i{x) < —1 for —1/e < X < 0. In addition, we require 

that Wn{-^) < —1 in order to have k(z; 13) > 0 in (191b . 
a^e 

Therefore, n = —1 is the only acceptable branch that results in 
real positive values of k(z; (3). This provides another proof for 

the fact that —1/e < < 0 or equivalently —a^ < /3 < 0 

_ cr^e 

(see Lemma 12J). 

From dglT l. for every fixed value of [3, we have 

= k(z;^) X ilog2 ^ 1 -f , z>0. (92) 

On the other hand, this choice of /? results in the instantaneous 
distortion D{z) := D{p{z), k{z)) in (|7]i as follows 

\ /3) X i log 2 (^1 + , z > 0. (93) 

Therefore, equating the l.h.s of ( l92l and ( l93b yields 



z > 0. (94) 
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Let /3* denote the optimal value of the free parameter jS. Thus, 
( l94l i results in the following associated optimal instantaneous 
distortion D*{z) 

D*(U = I “P (^) +0 * ^ “ (95) 

[ 0-2 z = 0 . 

As elaborated in Section I1V-B2I and Appendix |B] ( |95] ) is the 
unique constant solution of (l76l l for ^ > 0 in the Gaussian 
case. Specifically, p[z) and k{z) are jointly optimized such 
that when using a high (respectively low) transmission power, 
the transmitter uses a low (respectively high) mismatch factor 
to maintain the optimal instantaneous distortion to a constant 
value. 


VI. Numerical Results 


In this section, we consider numerical solutions to dhSl l and 
(|76]l to obtain an efficient power policy p{z) and mismatch 
factor policy k{z). In particular, although there is no apparent 
closed-form solution to dMI) for p{z), for every choice of the 
constants Gi, C 2 and the initial condition p(0+), we can apply 
numerical ODE solution methods. More precisely, given a fixed 
choice of (3, from ( f76b we can in principle solve for k.{z) in 
terms of p{z) (although other than the Gaussian case where a 
closed-form is found in dMT ), this must be done numerically). 
We then substitute k{z} thus computed into (l68l l and obtain an 
ODE for p{z) in terms of Ci, C 2 and /3 that can be solved 
numerically. Once p{z) is thus found, one can then directly 
obtain k{z) using (|76l) again. 

We obtain from (l20b that 


_ S 

TTo -f /o+ e^'^f{u)du p{z) + £{z) ’ 


(96) 


where by integrating both sides over (O+jZ] and performing 
some simplifications, this becomes 


Tio + / = TToexp [ [ 

Jo+ \Jo^ 


f , du . (97) 

lo+p{u)+i{u) J 

Taking the derivative of both sides in Wt\ with respect to z, 
we compute the density f{z), provided that ttq is known, as 
follows 


f{z) = TTo 


Se 


— Xz 


p{z)+i{z) 


exp 


' 0 + piu)+iiu) 


du 


(98) 


By combining ( l98T l and (l32l i. we compute the atom ttq as below 

\ -1 


TTo = 1 


(5e 


— Xz 


I 0 + piz)+iiz) 


exp 


I0+ p{,u)+i{u) 


du dz 


(99) 


Moreover, the value of the mismatch factor k(0) := k(z)|z=o 
when the battery is exhausted is obtained from ([33]) as 


k ( 0 ) =- j - — -. ( 100 ) 

1 -/ 0 + f{z)lKiz;l3)dz 

We now study a single-user EH communication system with 
energy arrival rates <5 = A = 1, and noise power N = 1. To find 
locally optimal policies for p{z) and k{z), one needs to search 
for optimized values of the free constants /3, Ci, and G 2 . In 
the following we separately investigate the cases of Gaussian 
and binary sources. 


A. A Gaussian Source over a Gaussian Channel 

Eor a standard Gaussian source 1), we have the bound 
— 1 < /? < 0 from Section IIV-B2I We now examine two 
cases of leakage; (i) zero leakage rate t{z) = 0 for an ideal 
battery and {ii) increasing leakage rate £{z) = 1 — for an 
imperfect battery. Obviously, our analysis is general enough to 
allow us to study different leakage behaviours i{z). However as 
stated in ED, “batteries leak most right after being charged", 
and examples of rechargeable batteries with leakage (or self¬ 
discharge) rates that increase monotonically with the battery 
charge are nickel-cadmium and nickel-hydrogen cells E9l . 

We also consider the battery capacities L = 1,2,..., 5. Table 
U and Table |II| show the total average distortion Davg, the 
distortion lower bound I?lb, and good values of the constants 
/3, Cl, and C 2 found by numerical search for both cases. 
Eurthermore, the initial condition of the ODE for p{z) was 
chosen to be = 0.001. This choice of p(0“'') is justified 

by the fact that a small amount of available energy in the 
battery should entail a small transmission power, as otherwise, 
the battery will be completely depleted before the next energy 
arrival. Numerical simulations have shown that if p(0’'") is 
chosen to be sufficiently small, it does not significantly change 
the distortion performance. A similar observation was made in 

El- 

It is evident from Table HI that as the battery capacity 
increases, the achieved distortion decreases. In particular, for 
L = 5 when the leakage rate is zero this scheme can achieve 
a distortion that is at most 7% above the lower bound. In 
fact, for an ideal battery with infinite capacity as discussed 
below ([38]), the lower bound is asymptotically tight and can 
be approximated arbitrarily well with a constant transmission 
power policy and a constant mismatch factor policy. Moreover, 
for case {ii) with non-zero leakage, the achieved average 
distortions are larger compared to that of case (t). Note that 
the distortion lower bound in Table HU is relatively more loose, 
since it does not depend on the leakage rate (see (l37l) ). while 
for the case of a non-zero leakage the optimal performance 
should depend on the leakage model. Therefore, we do not 
expect the lower bound to be asymptotically achievable, as 
it is without leakage. As the leakage model has increasing 
leakage with increasing battery charge, it is not unreasonable 
for the lower bound to become more loose with larger battery 
capacities. Nevertheless, the ratio in Table HJ must saturate to a 
finite value as L — c ». This is because i) the distortion lower 
bound for an infinite battery capacity is ((5/A, 1) which is 
strictly positive in this case, and ii) the achieved distortion will 
converge to some value, and hence, so will the ratio. 

Fig- ID shows the transmission power policy p{z), for the 
case of L = 5, for an ideal battery as well as an imperfect 
battery with three different leakage rates, i.e., increasing £{z) = 
1 — e“^, decreasing £{z) = e~^, and constant £{z) = 1. 

We observe that in all cases, the designed transmission 
power monotonically increases as the battery charge increases. 
This is due to the fact that when the remaining charge in the 
battery is close to the capacity limit, new energy arrivals are 
likely to make the battery overflow. Therefore, the transmitter 
consumes a large transmission power in order to avoid losing 
energy. Interestingly, for the increasing leakage and the con- 
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TABLE I: Distortion lower bound Dlb. average distortion D^vg, and good values of the constants Ci, C 2 and /3 of a standard Gaussian 
source A/'(0,1), for different battery capacities with the initial condition p(0’'") = 0.001, when (.{z) = 0. The ratio of Davg to Dlb quantifies 
the gap between the average distortion and the lower bound. 


Capacity of the Battery 

Dlb 

Davg 

-C^avg/ -f^LB 


Constants 


L = f 

0.6127 

0.6971 

1.13 

fi = 

-0.9485, Cl = -0.95, C2 

= 0.24 

L = 2 

0.5363 

0.6147 

1.14 

li = 

-0.9137, Cl = -0.92, C2 

= 0.30 

L = 3 

0.5128 

0.5765 

1.12 

fi = 

-0.8940, Cl = -0.90, C2 

= 0.32 

L = 4 

0.5046 

0.5559 

1.10 

fi = 

-0.8822, Cl = -0.89, C2 

= 0.32 

L = 5 

0.5017 

0.5417 

1.07 

fi = 

-0.8738, Cl = -0.89, C2 

= 0.34 


TABLE II: Distortion lower bound Dlb, average distortion Davg, and good values of the constants Ci, C 2 and /3 of a standard Gaussian 
source 1), for different battery capacities with the initial condition p(0’^) = 0.001, when l{z) = 1 — . The ratio of Davg to Dlb 

quantihes the gap between the average distortion and the lower bound. 


Capacity of the Battery 

Dlb 

Davg 

Davg/ 


Constants 


L = 1 

0.6127 

0.7445 

1.21 


-0.9641, Cl = -0.97, C2 

= 0.06 

L = 2 

0.5363 

0.6876 

1.28 

fi = 

-0.9450, Cl = -0.95, C2 

= 0.17 

L = 3 

0.5128 

0.6659 

1.29 

fi = 

-0.9366, Cl = -0.94, C2 

= 0.29 

L = 4 

0.5046 

0.6596 

1.30 

fi = 

-0.9300, Cl = -0.93, C2 

= 0.32 

L = 5 

0.5017 

0.6566 

1.30 

fi = 

-0.9302, Cl = -0.93, C2 

= 0.34 



Eig. 2: Power policy p{z) under 4 different scenarios, namely, zero 
leakage, increasing and decreasing leakage and constant leakage, 
when p(0''') = 0.001. 


slant leakage cases, the allocated transmission power increases 
faster compared to the transmission power for an ideal battery 
with the same battery charge. This result is intuitive, since 
an efficient transmission power policy mitigates the potentially 
large energy loss due to leakage by rapidly consuming the 
stored energy before it is lost. Eigs. [3 and |4] illustrate the 
corresponding mismatch factor k{z) and absolutely continuous 
part of the density function of the available charge in the 
battery, respectively. It can be seen that as the energy in 
the battery decreases, the mismatch factor increases. In other 
words, the low transmission power due to reduced charge in 
the battery is compensated by using longer channel codewords. 
Conversely, when the transmission power is large, codewords 
of smaller length are used so that the constraint in (l24l) is 
satisfied. As expected, an increasing transmission power results 
in a decreasing density function f{z). In other words, the 
storage process spends a smaller fraction of time at higher 
battery charges that have larger transmission power. 



Eig. 3: Bandwidth mismatch factor k{z) under 4 different scenarios, 
namely, zero leakage, increasing and decreasing leakage and constant 
leakage, when p(O^) = 0.001. 


The average distortion for the two cases of dynamically 
varying bandwidth mismatch factor and constant bandwidth 
mismatch factor is illustrated in Eig |5] both as functions of 
the battery capacity. It can be observed that for a battery with 
capacity in the range 2 < L < 12, a communication system 
with an adaptive mismatch factor, as proposed in this paper, 
performs better compared to that of k{z) = 1. Whether this 
gap is considered small or not depends on the application. 
If it is negligible for some applications, then using a fixed 
channel coding scheme would result in a reduced complexity 
system with negligible impact on performance. On the other 
hand, regardless of whether we are using a dynamic mismatch 
factor or a constant mismatch factor, for large values of L 
the average distortion of both coding schemes approach the 
lower bound and merge asymptotically. Likewise, for small 
battery capacities the average distortion of both coding schemes 
approach Umax and merge as well. 
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TABLE III; Distortion lower bound Dlb, average distortion Davg, and good values of the constants Ci, C 2 and /? of a Bernoulli(l/2) 
source, for different battery capacities with the initial condition p(0’^) = 0.001, when (.{z) = 0. The ratio of Davg to Dlb quantifies the gap 
between the average distortion and the lower bound. 


Capacity of the Battery 

Dlb 

f^avg 

f^avg/ 


Constants 


L = 1 

0.1651 

0.2097 

1.27 

/f = 

-0.3450, Cl = -0.36, C 2 

= 0.13 

L = 2 

0.1270 

0.1663 

1.30 

l^ = 

-0.3170, Cl = -0.32, C 2 

= 0.15 

L = 3 

0.1161 

0.1473 

1.26 

/f = 

-0.3039, Cl = -0.31, C 2 

= 0.16 

L = 4 

0.1122 

0.1367 

1.21 

/f = 

-0.2962, Cl = -0.30, C 2 

= 0.16 

L = 5 

0.1108 

0.1321 

1.19 

fi = 

-0.2901, Cl = -0.29, C 2 

= 0.16 


TABLE IV: Distortion lower bound Dlb, average distortion Davg, and good values of the constants Ci, C 2 and /3 of a Bernoulli(1/2) 
source, for different battery capacities with the initial condition p(O^) = 0.001, when l{z) = 1 — e~^. The ratio of Davg to Dlb quantifies 
the gap between the average distortion and the lower bound. 


Capacity of the Battery 

Dlb 

Davg 

Davg/ 


Constants 


L = 1 

0.1651 

0.2356 

1.42 

/f = 

-0.3600, Cl = -0.36, C2 

= 0.04 

L = 2 

0.1270 

0.2044 

1.60 

fi = 

-0.3417, Cl = -0.35, C2 

= 0.11 

L = 3 

0.1161 

0.1935 

1.66 

fi = 

-0.3347, Cl = -0.34, C2 

= 0.15 

L = 4 

0.1122 

0.1905 

1.69 

fi = 

-0.3328, Cl = -0.34, C2 

= 0.17 

L = 5 

0.1108 

0.1885 

1.70 

fi = 

-0.3301, Cl = -0.33, C2 

= 0.18 



- Dynamic mismatch factor 

- ' Constant mismatch factor k(z)=1 

— Lower bound 


0.8 -i 

I 

t 
% 
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10 15 20 
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25 


30 


Eig. 4: Density function f(z) under 4 different scenarios, namely, 
zero leakage, increasing and decreasing leakage and constant leakage, 
when p(0''') = 0.001. 


Eig. 5: Average distortion for the two cases of adaptive mismatch 
factor and constant mismatch factor k(z) — 1, when the battery 
capacity varies in the range 0 < L < 30. 


B. A Binary Source over a Gaussian Channel 


We now consider a binary source with the Bernoulli(1/2) 
distribution for which we have Lf(p) = 1, and the bound 
—0.5 < /3 < 0 as explained in subsection IIV-B2I Analo¬ 
gous to the Gaussian source, with different battery capacities 
L = 1,2,..., 5 and for the two cases of (i) zero leakage rate 
^{z) = 0 and {ii) increasing leakage rate i{z) = 1 — e~^, we 
evaluate the achieved average distortion Davg, the distortion 
lower bound Dlb and good values of the constants found by 
numerical search. The results of case (*) and {ii) which are 
summarized in Tables [HI] and respectively, have the same 
trends as in Tables HI and HIl with increasing the battery capacity. 
Moreover, the numerical results further show that for L = 5 
the proposed achievable scheme with no leakage can achieve 
a distortion which is at most 19% above the lower bound. 


VIE Conclusion 

We have investigated the problem of joint source-channel 
coding in a point-to-point channel with an energy harvesting 
transmitter. We used a calculus of variations technique to 
characterize an achievable joint source-channel coding scheme 
as well as an achievable transmission power policy to minimize 
the distortion at the receiver. We also obtained a distortion 
lower bound, where we used the convexity of the distortion 
function and an upper bound on the average transmission 
power. 

Eor a moderate-size battery capacity, we numerically showed 
that the achievable distortion with a dynamically varying 
bandwidth mismatch factor is smaller than that of a constant 
mismatch factor when the battery has no leakage. Moreover, 
we observed that as the battery capacity tends to infinity the 
achievable distortion for both coding schemes approached the 
lower bound. Eurthermore, with a constant mismatch factor 
k{z) = 1, we found a constant transmission power policy that 
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can arbitrarily approach this lower bound for an infinite battery 
capacity. 

As examples of continuous and discrete alphabet sources, we 
considered both Gaussian and binary sources to validate our 
analytical findings. In both cases, we showed numerically that a 
good transmission power policy increases as the battery charge 
increases. In contrast, a good mismatch factor policy, which 
measures the ratio of the length of channel codewords per 
source symbol, is a decreasing function of the battery charge. 
We further examined these policies under different possibilities 
for battery leakage rate, i.e., zero leakage rate as well as non¬ 
zero arbitrary leakage rate. 

One possible future extension is to consider a scenario where 
lossy source-channel communication is carried out over a time- 
varying channel with slow fading. In particular, the transmis¬ 
sion power and mismatch factor in this case need to be adapted 
to both the battery charge and the channel state. However, 
similar to Eol, we expect a system of coupled ODEs with one 
ODE per channel state, which can be challenging to solve, even 
numerically. Other future work includes considering other and 
more general energy arrival models. 


Appendix A 

To prove the convexity of the distortion function {p, q) 
over the domain p > 0,q > 0 where q := 1 /k, we consider two 
cases; (a) Rf' = oo, (b) Rf' < oo. As already discussed, an 
example of the former is a Gaussian source and an example of 
the latter is a binary source. We first prove the joint convexity 
with respect to p and q for case (a). To do so, we compute the 
Hessian matrix of D^{p^ q) for p > 0, g > 0, denoted by Hd, 
and show that it is positive definite. The Hessian is given by 


d^D^\ 



\ dq dp dq^ j 

where by a simple calculation we obtain 




dp'^ 


(a) 

> 


K{p) 

R'AD{p,k)) 


+ iKip)f X 


-R's {Djp, A) 

mDip,K))f 


0 , 


where (a) follows from the conditions [S2] and [C2] on 
Rs{D) and Rc{p), respectively. Moreover, the strict inequality 
is justified by the fact that since R}^ = oo, the distortion D 
is strictly positive and therefore R's{D) is finite due to [S3]. 
Similarly, we can show that the determinant of Hd is strictly 
positive, i.e., 

detiHo) = {Rc{p)f Ar'Ap) 

1 -RAD{P.A) A n 

R',{D{p,k)) RAD{p,A) 

where (b) again follows from conditions [S2], [S3] and [C2]. 
Since, iTuii > 0 and det{Hjj) > 0, by Sylvester’s criterion 
the matrix Hjj is positive definite, and it thus follows that 
Z)^(p, q) is jointly convex over the pair p and q. 

To prove the joint convexity with respect to p and q for 
case (b) where Rf^ < oo, we recall that {p, q) = 0 for 


Rc{p) > qRf'- Although the function D\p,q) is continuous 
everywhere, and in particular at the points where kRc{p) = 
RA, the second derivative at these points may not necessarily 
exist and a more complicated analysis is required. Therefore, as 
illustrated in Eig. |6]we separate the region p > 0 and q > 0 into 
two parts: the open shaded region, Rc{p) < qRf', over which 
the convexity argument reduces to the case (a) and the closed 
unshaded region, Rc{p) > qR^A which D^{p,q) = 0. It 
is not hard to see that the unshaded region given by Rc{p) > 
qRf^ (or equivalently q < Rc{p)/Rf') is convex. This is due 
to the fact that Rc{p) is a concave function and thus the region 
it traces is convex. Now, consider two arbitrary points = 
{pi,qi) and a 2 = (^ 2 , 92 ) such that ai,a 2 € {{p,q) ■ P > 
0,q > 0}. Eor A G [0,1], we define the function g{X) as 
follows 

g{X) := 1?^ (Aofi + (1 - A)q; 2 ) 

= D\Xpi -b (1 - X)p2,Xqi -b (1 - X)q2). 

If g{X) is convex for all choices of ai and 0 : 2 , so is D'^{p, q). 
With respect to the closed line segment connecting these two 
points, i.e., C = {Aai -b (1 — A)a 2 : A G [0,1]}, the following 
different cases may happen: 

• If ai and 0:2 are both in the unshaded region, the line 
segment C only passes through the unshaded region as 
the region is convex and thus g{X) = 0, which is convex. 

• If ai and 0:2 are in two different regions (say ai is in 
the unshaded region and a 2 is in the shaded region), there 
exists Ai G (0,1] as shown in Eig. |2h such that £ lies in 
the shaded region for A < Ai and it enters the unshaded 
region for A > Ai. Eurthermore, once the line segment 
enters the convex unshaded region, it does not exit. Here, 
g(X) is continuous for A G [0,1], non-negative and strictly 
convex for A G [0, Ai), while g(X) = 0 for A G [Ai, 1]. 
Thus, g(X) is convex. 

• If ai and a 2 are both in the shaded region, then either £ 
only passes through the shaded region where the convexity 
of g(X) reduces to the case (a), or due to the convexity 
of the unshaded region it enters the unshaded region for 
one and only one contiguous closed interval [A 2 , A 3 ] C 
(0,1) and again returns to the shaded region for A > A 3 . 
The function g(X) in the latter case is again continuous 
for A G [0,1], non-negative and strictly convex for A G 
[0, A 2 ) U (A 3 ,1], while g(X) = 0 for A G [A 2 , A 3 ]. Thus, 
g(X) is convex and this case is illustrated in Eig. [TJi. 

• Similar to the second case, when ai is in the shaded 
region and 0:2 is in the unshaded region, we have Eig. 

where g(X) is convex. 

Consequently, even for sources for which zero distortion can 
be achieved (e.g., the binary source), the function D'^{p,q) is 
jointly convex over the pair p and q. 

Appendix B 

We first rewrite (|7]i for z > 0 as 

Rs{D{p{z),k(z))) = k{z)Rc{p{z)), 


( 101 ) 
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Fig. 6: The shaded region shows the area where q > Rc{p)/Rf". 





Fig. 7: Special cases of the function g{\). 


where by taking the first derivative of both sides with respect 
to k{z) at a fixed z > 0 we obtain 


dD{p{z), p {p{z),k{z))) = Rc{p{z)), (102) 


dK{z) 

or equivalently 

dD{p{z),K{z)) 


Rc{p{z)) 


dK{z) R',{D{p{z),k{z)))' 


(103) 


To simplify the notation, we fix z > 0, and simply write D in 
place of D{p{z), k{z)). Thus, with the substitution (I1031 l. (l76l l 
reduces to 


D- 


RsjD) 
R's (D) 


+ /3 — 0, 


(104) 


where for the second term we also used the substitution (IIOII) . 
Therefore, D must be a root of (11041) . Solving (11041) for f3, we 
then obtain 


Rs{D) 

R's (D) 


-D. 


(105) 


We next show that the r.h.s of (11051) is strictly decreasing with 
respect to D. To do so, we take the first derivative of the r.h.s 
in (1105b with respect to D, and show that it is always negative 
in the open interval 0 < D < i.e.. 


d \Rs{D) 

dD [R', {D) 


-R'^{D)Rs{D) (c) ^ 

{R's{D)f 


(106) 


where (c) follows from the assumption of R'l{D) > 0 and 
further recalling that Rs{D) > 0 and R's{D) is finite, for 0 < 
D < Umax- Therefore, for every fixed /3, there is at most one 
real root D that solves (11051) . and it does not depend on z. This 
is also true for the optimized value of /?*, and therefore the 
associated instantaneous distortion D* {p{z), k{z)) is constant 
for z > 0. Furthermore, for there to be at least one real root, 
/3 in (11051) must be in the range /3niin < (d < /3max, where 


sup 

0<£l<D„,ax 


= lim 

DiO 


R. 


Rs{D) 

[R'siD) 

{D) 


-D 


R'sioy 


(107) 

(108) 


and 


inf 

<Z)niax 

lim 

Dy Dmax. 


Rs{D) 

[R'siD) 

Rs{D) 

R'siD) 


- D 


-D„ 


(109) 

( 110 ) 
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